Collective Suppression of Linewidths in Circuit QED 
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We report the experimental observation, and a theoretical explanation, of collective suppression 
of linewidths for multiple superconducting qubits coupled to a good cavity. This demonstrates how 
strong qubit-cavity coupling can significantly modify the dephasing and dissipation processes that 
might be expected for individual qubits, and can potentially improve coherence times in many-body 
circuit QED. 



At the root of many of the unexpected effects predicted 
by quantum mechanics is quantum interference. In the 
context of quantum optics, one notable consequence of 
constructive interference is superradiance 0, Q , the col- 
lective enhancement of radiation from an ensemble of 
many initially excited atoms. At its heart is the idea 
that if the atoms remain in a symmetric state, then con- 
structive interference between different final states after 
emitting a photon can increase the probability of such 
photon emission events It is notable that such con- 
structive interference plays a role even when considering 
the incoherent and irreversible emission of photons into 
an external environment. This collective effect is distinct 
from the subnatural linewidth averaging first discussed 
by 1. 

Since the superradiant emission of photons can oc- 
cur incoherently, it can be expected that similar ef- 
fects should be visible in the linewidth of a collection of 
atoms, or artificial atoms, coupled symmetrically to the 
same solid state environment. Indeed, one may see such 
an effect by comparing collective and individual decay 
processes. Consider N two-level systems (2LS) obey- 
ing either individual decay and dephasing: p — pn + 
Ei(l\\/2Maf] + [ lx /2)V[cxr] (where V[X] = 2XpX^ - 
X'Xp — pX^X, {<Ti} are Pauli operators describing the 
two-level systems with associated dephasing (711) and re- 
laxation (7j_) rates, and pn = —i[H,p] describes the 
Hamiltonian evolution), or collective decay and dephas- 

ing: p = p H + ( 7 ||/2)I?E l of] + (7±/2)2?E j H one 
then calculates the linear absorption spectrum for an en- 
vironment coupling symmetrically to all 2LS one finds a 
total linewidth I/T2 = 2711 + 7j_/2 for the case of individ- 
ual decay, and 1/T 2 = 2711 + Nj±/2 for collective decay. 
The linewidth associated with the coupling to a common 
bath is collectively enhanced, because of the constructive 
interference of different decay pathways. 

A natural context in which such questions arise is solid 
state realizations of coupled matter-light systems, such 
as circuit-QED 0, [Bj], where multiple superconducting 



qubits can be confined in a single microwave cavity 
and so may potentially couple to a common reservoir. In 
the limit of a good cavity, where a significant part of the 
vacuum Rabi linewidths is due to non cavity-mcdiatcd 
decay and dissipation, the distinction between coupling 
to collective and separate decay channels should be ap- 
parent in the dependence of linewidth on the number of 
qubits present. Even within a single sample, the effec- 
tive qubit number can be easily varied by detuning the 
qubits away from resonance with each other jfj. This 
breaks the symmetry, providing which-path information, 
and thus destroys the coherence. The question of whether 
decay and dephasing of multiple qubits is due to sepa- 
rate or collective coupling to the environment may have 
important consequences for the ability to preserve and 
manipulate coherence. For example, unexpectedly long 
coherence in light-harvesting-complexes 0] is associated 
with non-trivial quantum d yna mics arising from coupling 
to common photon modes [10| . 

There is however a problem with the simple picture of 
collective enhancement of linewidth when applied to mul- 
tiple qubits coupled to a microwave cavity. The problem 
is that the Lindblad terms written in the above are those 
that would be derived by considering system-reservoir 
coupling where the system Hamiltonian is that of a sin- 
gle qubit. The importance of using the correct system 
Hamiltonian in deriving loss terms has long been recog- 
nised in the context of ensuring that the correct equilib- 
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rium state is reached asymptotically [ll|, [12| . More re- 
cently there has been significant activity on such issues in 
the context of quantum dots [l3l - l22j | . In particular, it has 
been noted that when the system Hamiltonian is signifi- 
cantly changed, either by external driving |14l422| . or (as 
in the current case) strong qubit-cavity coupling|l3l. fl9j | , 
it is crucial to recalculate the decay processes in the pres- 
ence of the full system Hamiltonian. This is particularly 
important when the reservoir has a non-flat frequency de- 
pendence. While a sufficiently rapidly varying frequency 
dependence may prevent a Markovian density matrix be- 
ing used at all, there is a significant range of parame- 
ters where a Born-Markov approach remains valid, but 
it is necessary to calculate the decay rates using the cor- 
rect system Hamiltonian, rather than regarding the decay 
rates as fixed parameters. It is clear that the system of 
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many qubits coupled to a common photon mode is such 
a case: When the coupling and number of systems be- 
comes large enough, the system has been predicted Q 
to undergo a phase transition to a spontaneously polar- 
ized state 24|; at this point the frequency of the collec- 
tive mode should vanish, and care [25[ must be taken 
to avoid unphysical predictions. The importance of us- 
ing the eigenstates of strongly-coupled Hamiltonians has 
also been pointed out in the context of quantising super- 
conducting circuits 26]. 

In this letter we report measurements and calculations 
of the linewidth for one, two and three qubits coupled 
resonantly to a microwave cavity. In contrast to the 
well studied effect of superradiance we observe and ex- 
plain the narrowing of linewidths in the strong dephasing 
regime (i.e. dephasing and cavity decay rates are com- 
parable) as the number of qubits is increased. We find 
the results are compatible with a model of "collective" 
dephasing processes, where all qubits couple to a sin- 
gle bath, with a spectrum corresponding to 1// noise. 
We also discuss how varying the cavity-qubit detuning 
might allow corroboration of this scenario, and a direct 
measurement of the dephasing bath spectrum. 

Figure Q] shows part of the vacuum Rabi transmission 
spectrum (inset) and the extracted linewidths measured 
for a microwave cavity with one, two and three qubits 
tuned into resonance with the cavity mode. The super- 
conducting microchip sample and setup is similar to the 
one used in Refs. [a, [53]. The transmission spectrum 
is measured with much less than a single intra-cavity 
photon on average and clearly shows the expected \/N 
dependence of the vacuum Rabi splitting with increas- 
ing number of qubits (not shown). The single qubit- 
photon couplings of the three superconducting transmon 
qubits are g/(2n) = (52.7,55.4,55.8) MHz. The over- 
coupled coplanar waveguide resonator has a first har- 
monic resonance frequency of v r = 7.0235 GHz and 
a quality factor of Q = 14800 as measured when the 
qubits are far detuned (corresponding to a cavity decay 
rate k/2tt = 0.47 MHz). Time-resolved off-resonant T\ 
and T2 measurements of the three qubits confirm that 
T\ ^> T2 ~ 150 ns and that the qubit linewidths have 
strong dephasing components, i.e. I/T2 re for each 
qubit is fulfilled. The dressed linewidths on resonance 
are therefore expected to be dominated by qubit dephas- 
ing as well. While there is a considerable experimental 
uncertainty to the individual measurement data, see in- 
set in Fig.[T] the extracted linewidth from 6 single qubit, 
2 two qubit and 2 three qubit Rabi peaks (as indicated 
with error bars), shows a very clear trend of linewidth 
narrowing as the number of resonant qubits is increased. 

As discussed above, if the dephasing could be mod- 
eled by Lindblad operators derived for the uncoupled 
system, then narrowing of the linewidth would be un- 
expected. For a fixed dephasing rate, one would expect a 
constant or increasing linewidth. As discussed below, one 
can however directly explain this narrowing as a result of 
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FIG. 1. Experimental linewidth (black, solid) and theoretical 
linewidths (red and blue dashed) of an iV-qubit-cavity system 
(error bars indicate the sample standard deviation of the fitted 
linewidths). Theoretical linewidths are calculated for either 
collective coupling (lower line) or individual coupling (upper 
line) to a reservoir with a 1// spectrum (density of states 
J(y > 0) = 0.0105<?/^). The inset shows the experimental 
spectrum of the coherently scattered transmission amplitude 
for 1, 2 and 3 qubits, re-scaled and shifted so that peaks are 
centered at the same frequency for comparison. The lines are 
the corresponding Lorentzian fits. For every number N of 
qubits, both Rabi peaks were measured and used to calcuate 
a single linewidth and its uncertainty. The inset only shows 
one set of Rabi peaks. 



coupling the full system Hamiltonian to a frequency de- 
pendent bath; we assume in the following that dephasing 
is due to a bath with a 1// spectrum 28]. The fre- 
quency at which the coupled system samples the bath 
depends on the energy differences between system eigen- 
states. These depend on the collective Rabi frequency, 
which scales as yf~N for N qubits. Thus, the effective de- 
cay rate decreases with increasing number of qubits. This 
simple argument explains the essential origin of the re- 
sults seen experimentally, however several complications 
occur when one actually calculates the effective linewidth 
using the strongly-coupled qubit-resonator Hamiltonian. 

As discussed also above, two different scenarios of de- 
phasing can exist; individual coupling to reservoirs and 
collective coupling. When accounting for the qubit- 
resonator coupling, the dephasing induced linewidth de- 
pends differently on N in these two cases. In the collec- 
tive coupling case, the linewidth scales as 1/vN as the 
simple argument given above suggests: The dominant ef- 
fect is the sampling the reservoir at the collective Rabi 
frequency. Small corrections occur due to the matrix el- 
ements between system eigenstates arising from coupling 
to the bath and the small but non-zero photon decay 
rates. In the case of individual coupling, 1/y/N is still 
the dominant effect, but as system eigenstates are now 
delocalized, a competing effect arises: The number of 
possible decay channels increases with qubit number, as 
cross-qubit terms induced by the resonator couple qubits 
to other baths. This gives rise to an initial increase, fol- 
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FIG. 2. Relative theoretical linewidth of a single qubit- 
cavity system as a function of resonator-qubit detuning. The 
linewidths of both Rabi peaks (squares: upper peak, circles: 
lower peak) decrease because the photon part has a small de- 
cay rate and the qubit part samples the 1/f decay bath at 
frequencies a/A 2 + g 2 . 



lowed by a decrease of the linewidth. Figure [T] shows 
the results of the linewidth vs number of qubits in these 
two cases of individual and collective decay, assuming a 
1/f spectrum for the reservoir. Further details of the 
calculation are given in the following. 

In order to disentangle the effects of the reservoir den- 
sity of states from the effects of the cross-coupling ma- 
trix elements, one may instead explore the dependence 
of linewidth on cavity-qubit detuning. For a single qubit 
coupled strongly to a cavity, the normal modes are super- 
positions of photon and qubit states. As the qubit-cavity 
detuning A is varied two effects occur: The nature of the 

modes changes, and the frequency at which the reservoir 

^2 xhe 



is sampled varies as the Rabi frequency a/A^ 
changing nature of the modes means that one mode be- 
comes more photon-like, and the other will have a large 
qubit weight. The changing Rabi frequency means that 
the reservoir is sampled at higher frequencies, and so for 
a 1// spectrum, the linewidth will decrease. The calcu- 
lated linewidth vs detuning for a single qubit is shown in 
Fig. [5] Both the changing nature of the modes (causing 
the linewidths of the two modes to differ) and the re- 
duction of the qubit dephasing occur on a similar scale, 
A ~ <?, as is clear from the figure. For multiple qubits, 
the detuning dependence of the linewidths still depends 
on the nature of the coupling (individual vs collective) as 
well as the density of states of the reservoir. Nonetheless, 
the detuning dependence of linewidth can corroborate a 
given model of dephasing, and provide clear information 
about the real spectrum of the environment which in- 
duces dephasing. 

We now turn to discuss in more detail how the dephas- 
ing rate is calculated using the Born-Markov approxima- 
tion [29j], with the full Hamiltonian. Using Pauli opera- 
tors <Ji' y,z to represent the qubits and bosonic operators 
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bi q are the bosonic modes of the environment with en- 
ergy f3i q whose coupling strength to the system is j q . 
The important quantity for the behaviour of the qubits 
is the combination Ji{v) — J2qliqd( l/ ~ The 



total Hamiltonian H 



-Hbath describes the 



strongly coupled qubit-resonator system and the dephas- 
ing term in the case of coupling to individual reser- 
voirs. For collective coupling one instead has -ffbath = 
(Ei O E g lq ( b \ + b q) + J2q Pq h q h q- In thc interaction 
picture the equation of motion is 
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P = YlJ dvJ ^ J di'{W')/^(<) - Pi{t)Pi(t')p] x 

1 — OO —OO 

f v (t - t') + [P l (t) P P l (t r ) - pPi^Piit)] f v (t' - i)| (3) 



with Pi(t) — e ays *ofe~ sya * for individual reservoirs 
and P(t) = e lffsyst erf) e _lffsyat for collective de- 
phasing. We assume the reservoir density of states to 
be J(v > 0) oc 1/v and f v (r) = /£(— t) = (n„ + 
Y)e~ %vr + n v e %vr depends on the Bose-Einstein occupa- 
tion ri„ which we may assume to be zero, since fc^T <C g. 

By assuming p = p(t) on the right hand side of Eq. 
([3]) we make the standard Markov approximation [2^], but 
retain the crucial dependence on the collective Rabi fre- 
quency by evaluating the remaining integrals with the full 
system Hamiltonian (l7j . For a general system Hamil- 
tonian, Eq. is not in Lindblad form and so does not 
necessarily preserve positivity of the density matrix equa- 
tion [2§| . For short times, this is not an issue 17 1, but 
for long times (as matters for steady states), positivity 
violation becomes a problem [2l], • One may however 
obtain a Lindblad form by dropping any terms which are 
time-dependent in the interaction picture 31|; the re- 
sulting time evolution is equivalent to Eq. (J3|) over short 
times 31|, but avoids unphysical positivity violations. 
The resulting Lindblad terms have the form: 



E ( 2 <sprt ~ r? a r> - pr* r* tf ) (4) 



with r % a o = A a fi\a) {a\af\f3) where a, (3 are eigen- 
states of the system Hamiltonian and A a p = r(e Q — ep) 
is a transition-energy dependent decay rate, T{8) = 
ir [J(S)(n(5) + 1) + J(—S)n(—5)]. The summation over 
states is restricted to energy-conserving transitions E a — 
Ep — Es — E 1 . For collective decay, the sum over i dis- 
appears, and r Q( 3 is defined in terms of matrix elements 
of the operator J2i °f instead. 
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This dephasing term is of Lindblad form as required 
to preserve positivity of the density matrix. However, 
it is not simply a collection of qubit dephasing terms 
add ed ad hoc to the density matrix equation of motion. 
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|New J. PhysTl 2, 083033 (2010)| A. W. Chin, J. Prior, 
R. Rosenbach, F. Caycedo-Soler, S. F. Huelga, and 
M. B. Plenio, |Nat. Phys. (2013) | IS SN 1745-2473. 



Inst ead, wo naturally obtain dephasing in terms of system 



[12 

eigenstates, sampling the reservoir density of states at 
the system's natural transition frequencies. To calculate j > 



the coherent photon scattering spectrum we numerically 



solve the steady state of this density matrix equation in 
the presence of a weak drive. (Since the drive is weak, 
it does not itself modify the decay rates). The spectrum 
is then found by calculating \(a)\ vs detuning, and the 
linewidth plotted in Figs. [T] and [2] is extracted by fitting 
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thislspectrum to a Lorentzian. 



In conclusion, we have shown how strong matter-light Fff 
coupling and a non-trivial reservoir spectrum can pro- 
duce a non-trivial suppression of linewidth, which can 
non etheless be explained within a Markovian approxi- __M 



mat Lob. — Our calculations explain why linewidths in an 
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resonant qubit-cavity system and offer a way to probe 
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the -eservoir density of states. 
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